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1.0  Introduction 

Increasing  demands  on  the  performance  of  satellite  communication 
antennas  have  generated  recent  interest  in  efficient  and  systematic 
methods  for  analyzing  and  synthesizing  reflector  antennas  and  associ- 
ated feed  systems.  A comprehensive  discussion  of  a diverse  array  of 
available  techniques  for  analyzing  reflector  antennas  has  recently  been 
published  by  Rusch  [1];  the  reader  interested  in  developing  a good  broad 
background  on  the  subject  will  do  well  to  refer  to  this  excellent  work. 

Scanning  through  the  literature  on  reflector  antennas,  one  finds 
that  most  workers  prefer  to  use  methods  based  on  integration  of  the 
induced  currents  on  the  reflector  to  compute  its  secondary  pattern. 

Also,  the  induced  current  on  the  surface  of  the  reflector  is  typically 
derived  via  the  well-known  "physical  optics"  approximation  which  has  been 
verified  to  be  extremely  accurate  for  predicting  the  pattern  over  a rather 
wide  angular  range  in  the  forward  direction. 

In  contrast,  the  ray  optical  methods  are  known  to  run  into  difficulty 
[2]  in  the  neighborhood  of  the  caustic,  i.e.,  the  beam  maximum,  although 
they  work  well  at  wide  angles.  Hybrid  techniques  that  employ  ray  methods 
as  a first  step  for  deriving  the  fields  in  the  projected  aperture  of  the 
reflector  and  then  integrating  over  this  aperture  to  derive  the  far  field 
have  been  used  by  a number  of  workers.  However,  they  suffer  from  the 
shortcoming  that  the  secondary  patterns  computed  in  this  manner  rapidly 
become  inaccurate  as  one  moves  away  from  the  beam  maximum.  Thus,  the  use 
of  this  approach  is  limited  to  the  computation  of  the  main  beam  and  the 
near-end  sidelobes  only.  Furthermore,  the  polarization  information  de- 
rived from  either  the  ray  or  the  hybrid  approach  can  contain  errors 
beyond  the  tolerance  requirements  on  the  accuracy  of  the  results. 

For  the  reasons  outlined  above,  we  choose  to  follow  the  induced 
current  integration  method  in  this  work  and  describe  in  detail  a recently 
developed  series  approach  [3]  that  allows  an  extremely  efficient  evalua- 
tion of  the  radiation  integral  of  reflector  antennas.  The  series 
approach  differs  both  in  spirit  as  well  as  in  detail  from  the  various 
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numerical  algorithms  conventionally  employed  for  increasing  the  efficiency 
of  the  numerical  procedures  [1]  for  evaluating  the  radiation  integral. 

The  development  of  the  series  representation  is  motivated  by  the  fact  that 
in  the  course  of  synthesizing  contour  beams,  a topic  which  we  discuss  in 
the  latter  part  of  this  paper,  one  finds  it  necessary  to  perform  accurate 
computation  of  the  secondary  pattern  of  reflector  antennas  for  a large 
number  of  laterally  displaced  feeds  and  over  a wide  range  of  observation 
angles.  The  series  representation  of  the  radiation  integral  developed  in 
this  paper  allows  one  to  accomplish  this  task  in  a much  more  efficient 
manner  than  would  be  possible  via  any  of  the  conventional  approaches. 

The  organization  of  this  paper  is  as  follows.  In  Section  2 we 
discuss  the  case  of  symmetrical  reflectors  in  some  detail.  In  Section  3 
we  deal  with  the  offset  reflector  and  show  that  certain  transformations 
incorporated  in  the  radiation  integral  for  this  case  allow  one  to  convert 
it  into  a form  which  is  identical  to  the  symmetric  case,  thus  preserving 
all  of  the  useful  properties  of  the  series  form.  In  Section  4 we 
discuss  the  problem  of  contour  beam  synthesis  using  a cluster  of  feeds. 

We  develop  two  approaches  for  deriving  the  excitation  coefficients  for 
the  feeds,  one  based  on  an  approximation  method,  and  the  other  requiring 
the  solution  of  a matrix  equation  that  can  be  inverted  in  a closed  form. 

A comparison  of  these  two  methods  is  also  included  in  Section  4. 

2.0  Evaluation  of  Radiation  Integral  for  Symmetric  Reflectors 
2. 1 Integral  Representation  of  Far  Field 

The  far  field  radiation  integral  (suitably  normalized)  representing 
the  field  produced  by  the  induced  current  distribution  j on  the  surface  of 
the  reflector  can  be  expressed  as 


E{e,<^)  = (i  - RM 


where  I = the  unit  dyadic. 


' J I 

SURFACE 


. -jk(p'-P-Rn) 

Je  ^ ds. 


(2.1) 


The  geometry  of  the  reflector  antenna  system  including  the  feed  is 
shown  in  Fig.  1.  We  point  out  that  the  optical  phase  of  j has  been 
factored  out  in  (2.1)  where  this  phase  is  defined  relative  to  an  ideal 
phase  center  at  e,  which  defines  the  location  of  the  displaced  feed 
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Fig.  1.  Feed -ref Tec tor  geometry. 

From  Fig.  1 we  also  have  the  relationship 


where  p is  the  distance  from  the  origin  (focal  point)  to  a point  on  the 
reflector. 


We  note  that  the  integral  (2.1)  is  defined  on  the  surface  of  the 
reflector.  Often  a ray-optical  "approximation"  is  used  to  derive  an 
alternate  form  for  the  radiation  integral  which  requires  the  Fourier 
transformation  of  the  electric  and  magnetic  fields  on  a projected  aper- 
ture of  the  reflector,  rather  than  on  the  reflector  surface  itself. 
However,  as  mentioned  earlier,  this  introduces  errors  in  the  secondary 
pattern  calculation  at  observation  angles  away  from  the  main  beam  direc- 
tion and  is  therefore  not  suitable  for  our  purpose.  We  can,  nevertheless, 
derive  a mathematically  exact  form  of  (2.1)  which  is  also  an  integration 
over  the  projected  aperture  coordinates.  For  instance,  we  can  write  (2.1) 


J 


equivalently  as 

(•21'fa  -jk(p'-p*Rr|) 

F(0,4i)  = ■f(r,(t)')  e ^ r dr  6^'  (2.3) 

Jo  Jo 

a = radius  of  projected  circular  aperture 

where  F(0,(|))  is  related  to  a scalar  component  of  E,  say  the  x-component, 
and  f(r,(i)')  is  expressed  in  terms  of  3 and  the  Jacobian  relating  the 
elemental  area  on  the  reflector  surface  into  one  on  the  projected  circular 
aperture.  Explicitly,  we  have 

f(r,,^')  = J^(r.,^')/l  + . (2.4) 

and  for  a parabola  of  focal  length  f, 

F(r,<|)')  = J^(r,(^') 

The  function  f is  sometimes  referred  to  as  the  "equivalent  aperture" 
distribution.  It  should  be  realized  that  this  is  somewhat  of  a misnomer 
since,  physically,  f is  not  the  distribution  on  the  projected  aperture 
(in  contrast  to  the  one  that  can  be  derived  approximately  via  the  ray 
technique).  Since  (2.3)  is  exact,  if  f were  indeed  the  true  aperture 
distribution,  the  kernel  of  (2.3)  would  have  had  the  Fourier  transform 
form  which,  at  the  moment,  it  does  not.  One  of  the  principal  contributions 
of  this  work  would  be  to  show  that  the  kernel  in  (2.3)  can  still  be 
expressed  in  terms  of  a series  of  Fourier  transforms  of  appropriate 
functions  related  to  f. 

Typically,  the  induced  current  3 is  derived  from  the  physical 
optics  approximation 

3 = 2n  X H| 

where  H^.  is  the  primary  magnetic  field  incident  on  the  reflector  surface 
due  to  the  given  feed.  However,  if  desired,  the  physical  optics  approx- 
imation for  3 can  be  refined  by  augmenting  it  with  fringe  currents  existing 
at  the  edges  of  the  reflector. 
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Returning  now  to  (2.3),  we  note  that  the  conventional  approach  to 
evaluating  the  radiation  integral  entails  repeated  computation  of  this 
two-dimensional  integral  on  the  "equivalent  aperture"  for  each  observation 
direction  which  themselves  form  a two-dimensional  grid.  Even  if  some 
timesaving  can  be  arhieved  by  employing  efficient  numerical  algorithms  for 
evaluating  (2.3),  the  total  computation  time  can  become  prohibitively 
large  if  accurate  results  are  desired  over  a wide  range  of  observation 
angles.  This  fact  motivates  us  to  seek  alternate  means  for  handling  the 
radiation  integral  in  a manner  elaborated  on  in  the  following  section. 

2.2  Development  of  the  Series  Representation 

As  a preamble  to  developing  the  series  representation,  we  introduce 
certain  notations,  approximations  and  new  variables: 

p'  = p|l  - (^u-  +^v'  +^wjj  + 0(e2)  , (2.5a) 

where 

u’  = sin  9*  cos  ((>’;  v'  = sin  e'  sin  (f';  w'  = cos  e'  (2.5b) 

and  the  direction  cosines  u,  v,  w for  the  observation  variables  are  like- 
wise defined  by  removing  the  primes  from  (2.5b).  Typically,  for  small 
displacement  of  the  feed,  the  second-order  term  in  (2.5a)  is  negligible 
and  we  make  the  approximation  here  that  this  term  is  small.  However,  it 
j is  only  a trivial  extension  to  include  the  higher-order  term  in  (2.5a). 

Scaling  the  radial  coordinate  r by  introducing  r = as,  a = D/2,  and 
: D = diameter  of  the  reflector,  one  may  rewrite  (2.3)  as 

1 F(u,v)  ^ f'  [f(s.*') 

, Jo  ‘0 

■ . ^kass1necos(t-*')^  (2  5, 

I Upon  examining  the  integral  in  (2.6)  carefully,  we  find  that  the  first 

factor  inside  the  integral  is  a function  of  the  integration  variables 
alone  and  the  last  factor  is  the  Fourier  transform  kernel  written  in  a 
polar  form.  The  second  factor,  viz.,  exp{-jxpw' (1  - cos  0)},  is  a 
function  both  of  the  integration  and  observation  angles,  and  thus  is  the 
factor  that  requires  the  special  treatment  described  below. 
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We  define 


L(s,u,v)  = pw' (1  - cos  e) 


(2.7a) 


*-n  ~ L(s,Uq,Vq);  = L(Sg,u,v.;  = L(s^,Un»Vf,) 


aO 


(2.7b) 


where  is  a constant  to  be  determined  later  and  Uq,  Vq  are  the  direction 
cosines  evaluated  at  e = 0q,  (j*  = <)>q. 

Since  pw'  = p cos  0'  is  simply  the  z-coordinate  of  the  reflector 
surface,  for  a symmetric  parabola  this  quantity  is  a quadratic  function  of 
the  radial  variable  s,  and  is  obviously  indpendent  of  the  azimuthal 
variable  i(i'.  We  make  advantageous  use  of  this  property  of  L later.  If  we 
now  rewrite 


-jkpw*  (1-COS0)  _ ^-jL  _ 


= e ''  = e 


(2.8) 


we  may  rearrange  (2.6)  as 


E,(u,v)  = f'  f‘  . e-j‘.  s d*'  (2.9a) 

1 Jn  in 


where 


'0  ^0 


Ej  = F(u,v) 


(2.9b) 


f = fe 


(2.9c) 


a(u,v,Uq,Vq;  s.Sg)  = L - Lq  - Lg  + L^q  . 


(2.9d) 


One  last  step  in  our  manipulation  of  (2.6)  entails  the  shifting  of 
the  origin  of  the  observation  coordinates  such  that  it  is  centered  around 
the  beam  maximum  which  is  determined  by  I,  the  lateral  displacement  of  the 
feed.  This  transformation  is  useful  for  efficient  evaluation  of  the 
radiation  integral.  We  choose  the  new  observation  variables  to  be  n and 
a, related  to  the  original  direction  cosines  u and  v through  the  formulas 


7 2 ? * " ''M 

(u  - U|^)  + (v  - Vj,,)  ; a = arctan  ^ 

M 


and  being  the  coordinates  of  the  shifted  origin  which  is  here  taken 

to  be  the  beam  maximum.  Typically  Uq,  Vq,  introduced  in  (2.7b),  are  set 
equal  to  Uj,^,  Vj^,  respectively.  We  can  now  return  to  (2.9a)  and  rewrite 
the  first  exponent  as 


ke  • p + kas  sin  6 cos  (<()'  - 4i)  = ke,  cos  e*  + kas[-(B^^t„)  cos  (y|^  - 


'z '■"M'-M' 

+ n cos  (a  - 4)' )] 


where 


• I'M  " 


V V 

MM 


““M  ' '(V  * ') 


and  where,  for  a paraboloid, 

M H 1 + 


J_Y  s2 

D, 


(2.10) 


Here  we  have  made  use  of  the  relationshios 


~ 2(f/D)  M 


; cos  e = 


. - M - 2 


M • 


Lastly,  and  are  the  so-called  beam  deviation  factors  defined  by 


Bu  = 7^  (v  = 0 plane)  , B^,  = (u  = 0 plane) 


V V, 


with  u^  = . Simple  expressions  for  determining  B^^  and  B^  are 

given  in  Appendix  I.  Using  (2.10)  in  (2.9a)  we  can  rewrite  the  latter  as 


Ztt  1 

Ei(u,v)  = f f fi(s,({i’)  s ds  d<^'  (2.11a) 

^ Jn  Jn  ^ 


0 ■'0 


with 


fj(s,4>')  = f(s,(j)')  e e 


-jLf.  jke  cose'  -jkas6„t.,cos(Yw-().‘ ) 

^ e MM  M ^ (2.11b) 
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The  motivation  for  the  manipulations  leading  up  to  (2.11)  is  now 
given.  If  we  can  show  that  the  exponent  A is  small  for  certain  ranges 
of  observation  angles  to  the  extent  that  exp[-jA]  = 1,  we  would  obtain 
a Fourier  transform  form  for  the  radiation  integral  with  an  "effective 
aperture"  distribution  proportional  to  f^. 

Referring  to  the  definition  of  A given  in  (2.9d)  and  the  expressions 
for  the  various  L's  given  in  (2.7),  we  immediately  see  that  A = 0 at 
u = Uq,  i.e.,  0 = Sq  (typically,  though  not  necessarily,  the  elevation 
angle  corresponding  to  the  beam  maximum).  Thus  exp[-jA]  is  identi- 
cally equal  to  1 for  this  choice  of  observation  angles.  In  addition, 
though  not  obvious  at  this  point,  it  turns  out  that  by  a suitable  choice 
of  s,,  defined  in  connection  with  (2  7b),  we  can  make  A small  in  the 
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region  of  integration  of  (2.9a)  from  which  the  dominant  contribution  of 
the  radiation  integral  is  derived  in  the  wide  angle  region  (typically 
just  a few  si  delobes  away  from  the  beam  maximum)  where  e is  not  close  to 
0Q.  Thus,  unusual  as  it  may  appear  at  first  sight,  what  we  have  accom- 
plished through  the  manipulations  described  aoove  is  the  derivation  of 
an  "effective  aperture"  distribution  which  is  simultaneously  valid  not 
only  at  the  beam  maximum  (as  is  the  ray  optical  derivation),  but  at  wide 
angles  as  well  where  the  ray  optical  distribution  is  inaccurate.  In 
fact,  numerical  results  show  that  only  in  a small  intermediate  angular 
region  does  the  setting  of  exp[-jA]  = 1 in  (2.11)  introduce  any  signi- 
ficant errors. 

To  fill  in  this  gap,  we  go  now  to  the  next  step  and  expand  exp[-jA] 
in  a power  series  of  A,  obtaining 


2ir 


Ej(u,v)  = fjK  s ds  d(()'  = jj 


2w 


Af  jK  s ds  d(() ' 


r2iT 


0 Jo 


A fjK  s ds  d(t.' 


(2.12b) 


where 


K = gjkdSricos(a  (j>  ) _ transform  kernel  (2.12b) 


in  polar  form,  and  A was  defined  in  (2.9d).  At  this  point  we  can  compute 

2 

Ej  by  repeatedly  Fourier  transforming  f^,  Af^,  A f^,  etc.,  which  of  course 
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requires  an  n-fold  increase  in  the  number  of  computations,  where  n is  the 
number  of  terms  retained  in  the  series.  Instead,  in  the  next  section, 
we  introduce  a series  method  for  computing  the  first  integral  in  (2.12) 
which  has  the  unique  advantage  that  the  coefficients  of  the  series,  once 
computed,  can  also  be  used  to  evaluate  the  higher-order  terms  in 
(2.12)  with  very  little  additional  effort  or  increase  in  computational 
time. 

In  summary,  what  we  have  accomplished  in  expressing  the  radiation 
integral  in  the  form  of  (2.12)  is  a biconvergent  series  whose  leading  term 
accurately  predicts  the  pattern  both  in  the  neighborhood  of  e = Sq  and  at 
wide  angle  regions.  In  addition,  because  of  the  special  form  of  a,  the 
higher-order  terms  can  be  calculated  with  only  a little  additional  effort 
beyond  that  required  to  evaluate  the  first  term.  We  substantiate  these 
assertions  in  the  next  section  where  we  detail  the  evaluation  of  the 
various  Fourier  transform  integrations  in  (2.12)  using  the  Jacobi 
polynomials. 

2.3  Evaluation  of  (2.12)  Using  Jacobi  Polynomials 

The  first  integral  in  the  r.h.s.  of  (2.12),  which  has  strictly  a 
Fourier  transform  form,  can  be  evaluated  by  any  number  of  ways,  e.g., 
using  the  FFT  or  expanding  f^  in  a series  of  Bessel  functions  and  evalu- 
ating the  resulting  series  of  integrals  term-by-term  in  closed  form. 
However,  as  explained  in  the  last  section,  our  goal  is  to  utilize  the 
results  of  the  first  integration  for  efficiently  evaluating  the  higher- 
order  terms.  The  expansion  functions  for  f^^  that  are  uniquely  suited 
for  this  purpose  are  the  Jacobi  polynomials,  also  referred  to  as  the 
circle  polynomials  of  Zernicke  [4],  [5].  The  use  of  the  Jacobi  polynomial 
allows  one  to  evaluate  the  first  integral  with  a speed  comparable  to  the 
FFT.  However,  the  Jacobi  polynomial  method  has  several  important 
advantages  over  the  FFT  method,  viz.,  (i)  it  is  better  suited  for 
circular  apertures;  (ii)  it  can  be  used  to  compute  the  pattern  at  an 
arbitrary  value  of  e,  ())  not  restricted  by  the  choice  of  the  FFT-type 
grid;  (iii)  it  does  not  suffer  from  aliasing  errors  or  Gibb's  phenomenon; 
(iv)  it  is  useful  for  efficiently  evalu:  mg  higher-order  terms  whereas 
the  FFT  method  is  not;  (v)  it  has  superior  numerical  accuracy. 
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The  Jacobi  polynomial  series  expansion  begins  with  a two-dimensional  ' 

Fourier  series  representation  for  in  (2.12)  | 

j 

I 

A 

i 

(2.13)  • i 

i 

. i 

(2.14)  I 


fl(s,<f')  = I I fc  cos  n(^'  + dJJ  sin  n(^']  fJJ(s) 

^ m=0  n=0  ' 

f"(s)  s p"(1  - 2s^)  = 2(n  + 2m  + 1)  pi"’°^l  - 2s^)  s" 
m m m 


L 


where  pf,'^’*^^(x)  are  the  Jacobi  polynomials. 


Inserting  (2.13)  into  the  first  integral  in  (2.12),  henceforth 
referred  to  as  Ej  for  convenience,  and  using  the  identity 


2^  1 / \ ^ 
•IT  Jn(^)  = 

J 


2tt 

0 


cos  nu  du 


(2.15) 


we  get 


E^(n»a) 


„ , cos  na 

l2.jV  , 

n 1 sin  na  I m 


f"(s)  J^(kans)  s ds 


(2.16) 


where  denotes  an  upper  plus  a lower  sum.  The  integral  in  (2.16) 
involving  the  product  of  Jacobi  polynomials,  the  Bessel  function  and 
the  weight  function  can  be  evaluated  in  a closed  form  [6]  using  the 
integral  relationship 


1 ^nn/Z  p^n,0)(j  . 232,  ^ 


[n  > -(m  + 1)]  . 
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H*  III 


In  general  we  have 


I 


S ! 


'm.nl’')  = L = "s  = -^(n  * 2n,  * 1)  . (2.18) 


The  use  of  (2.18)  in  (2.16)  permits  us  to  write  the  final  form  for  the 
first  term  of  the  radiation  integral 


En^^(u,v)  = E(®^(n.a)  = 2it  j 
^ ^ n=0 


.n 


I 2(n  + 2m  + 1) 
m=0 


(2.19) 


(kan) 


.1 


We  note  that  the  radiation  integral  can  be  readily  calculated  at  any 
observation  angle  u,v  [or  (n.a)]  once  the  Fourier  coefficients  c[J|  and  0^^ 
have  been  determined.  This  is  in  contrast  to  the  direct  evaluation  of  the 
original  integral  (2.6)  which  has  to  be  computed  repeatedly  for  each 
observation  point.  We  also  note  from  (2.19)  that  the  leading  term  is  the 
Airy  function  which  is  the  radiation  pattern  of  a circular  aperture  with 
uniform  amplitude  and  phase  distribution  (except  for  a linear  taper  which 
shifts  the  beam  maximum  to  0^,1,  v^^).  All  of  the  higher-order  terms  of  the 
series  in  (2.19)  become  zero  as  n 0 and  represent  perturbations  around 
the  beam  maximum. 

The  expansion  coefficients  are  found  in  the  usual  manner  after 
recognizing  that  the  Jacobi  polynomials  satisfy  the  orthogonality 
relationships 


^ p"(l  - 2s^)  P".(l  - 2s^)  s ds  = 6"',  (2.20) 

Q m m m 


where  is  the  Kronecker  delta.  The  orthogonality  property  was  one  of 
the  criteria  employed  to  select  the  expansion  functions  for  f^;  the  other 
was  the  recursion  relationship  which  is  used  to  evaluate  the  higher-order 
terms.  The  integrals  to  be  computed  for  evaluating  the  series  coeffi- 
cients are  given  by 


fl(s,(li') 


cos  nij)' 
sin  nijj'i 


P"(l  - 2s^)  s d(j,'  ds 
m 


(2.21) 


where  = 1 fo>^  n = 0 and  2 otherwise.  An  important  characteristic  of 
the  integrands  in  (2.21)  is  that  they  are  far  less  oscillating  than  the 
kernel  of  the  radiation  integral  (2.6),  since  the  absence  of  the  exponen- 
tial factor  exp[jkas  sin  9 cos  ((p  - (j>')]  makes  the  kernel  of  the  radiation 
integral  in  (2.21)  vary  much  less  rapidly  compared  to  the  one  in  (2.6), 
particularly  when  {ka  sin  0}  is  not  small. 

Before  proceeding  with  the  discussion  of  the  steps  for  the  evaluation 
of  the  higher-order  terms,  we  examine  the  leading  term  Ej  in  a little 
more  detail  with  the  purpose  of  showing  that  this  term  is  useful  for  the 
wide  angle  region  as  well.  This  is  the  biconvergent  property  alluded  to 
earlier  in  connection  with  the  series  given  in  (2.12). 

In  order  to  explore  the  wide  angles  (e  not  close  to  e^)  behavior  of 
A,  we  refer  once  again  to  the  definition  given  in  (2.9d)  and  (2.7)  for  A 
and  the  L's,  respectively.  It  is  evident  that  A h 0 at  s = s . We  now 

d 

appeal  to  the  well-known  result  in  high-frequency  diffraction  phenomenon 
that  away  from  the  caustic  (beam  maximum)  the  contribution  to  the 
secondary  radiation  comes  mainly  from  the  rim  of  the  reflector,  described 
by  s = 1. 


Consequently,  if  we  choose  s^  to  be  equal  to  1,  we  make  the 
exponential  factor  e close  to  unity  in  the  neighborhood  of  the  region 
from  where  the  principal  contribution  of  the  radiation  integral  is 
derived,  in  the  asymptotic  sense,  for  all  observation  angles  not  close 
to  u = Ug  (e  = 0q).  Extensive  numerical  experiments  with  different 
choices  of  s,  have  confirmed  that  this  choice  for  s,  is  the  most 

a O 

appropriate  one  for  the  symmetric  reflector  antenna,  with  or  without 
displaced  feeds. 

The  last  topic  we  discuss  in  this  section  is  the  use  of  the  Jacobi 
polynomial  series  for  the  evaluation  of  the  higher-order  terms  in  (2.12). 
We  begin  by  explicitly  writing  L(s,u,v),  originally  defined  in  (2.7a). 

As  observed  earlier,  for  a parabola  pw'  is  simply  the  z-coordinate  of  the 
reflector  surface  and  this  variable  is  independent  of  if'  and  quadratic 
in  s.  Explicitly,  for  a paraboloid  of  revolution 


! Us.u.v)  = (.  ^)[.l  . (^)^]  sj}  . [.  a - f 

I . Using  the  definitions  given  in  (2.7b)  and  (2.9d),  we  can  write  A as 

t , A(L  - Lg)  = c s^  - Sg  G(ug,u)  (2.23a) 

f 

% 

where 
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Substituting  (2.23)  in  (2.12)  we  obtain 

Ej(u,v)  = I = e[°^(u,v)  + e[^^(u,v)  + Ep^(u.v)  + ...  (2.24) 

P 

where  E^^^,  defined  earlier, is  the  first  integral  in  (2.12),  and  the  higher- 
order  integrals  are  expressed  as 


E 


E 


(1) 

1 

(2) 

1 


s^j  flK  s ds  d({)' 


(2.25) 

(2.26) 


We  have  already  discussed  the  evaluation  of  E^®^  by  expanding  f^  in  a 
two-dimensional  Fourier  series  given  in  (2.13).  We  now  show  that  by 
virtue  of  the  recursion  relations  satisfied  by  the  Jacobi  polynomials, 
the  expansion  coefficients  of  the  higher-order  terms,  e.g.,  (s^  - s^J  f, 

(2  2\2  * n ” ' 

s - s^j  f^,  etc.,  can  be  readily  expressed  in  terms  of 
which  have  presumably  been  computed  already  in  connection  with  the 
evaluation  of  e|*^\ 

To  this  end  we  make  use  of  the  recursion  relationship  [6] 


where 


a P . + 
mn  m-1 


c P" 
mn  m+1 


m(n  + m) 

“mn  ~ ' (n  2m)(n  + 2m  + 1) 

u _ (m  n)^ (m  1)^ 

mn  ” (2m  + n)(n  + 2m  + 1)  (n  + 2m  + 2)(n+  2m  + 1) 


= _ (m  l)(n  1-  m+1) 
mn  (n  + 2m  + l)(n  + 2m  + 2) 


(2.27) 
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and  defining  the  higher-order  expansion  coefficients  in  terms  of  the 
integrals 


I 

I 


I 


'0  J 


Zt! 


e rl  ec-n  f 


COS  ni*' 


sin  n^i' 


p"(l  - 2s^)  s ds  (2.28) 


we  can  derive 


PrI  = 


m 


mn 


‘’■‘b"  , * 

m-1 


mn 


P-V  + 


s"  ■■  B ■ + c 
a m mn 


p-lnn 

Vl 


(2.29) 


where  B represents  either  C or  D.  We  note  from  (2.29)  that  if  we  retain 
terms  up  to  the  p^^  order,  then  the  number  of  coefficients  requiring 


computation  is  M + p where  M is  the  upper  limit  of  m.  Consequently, 


2N  (N  = n„,  ) additional  series  coefficients  (N  each  of  C's  and  D's)  are 
needed  for  evaluating  p integrals  of  the  type  Ej^'(u,v).  However,  one 
finds  that  the  higher-order  coefficients  typically  decrease  very 


rapidly,  allowing  truncation  to  a moderate  number  of  terms  and  correspond- 
ing pairing  of  computational  time. 


1 i 


We  summarize  this  section  by  noting  that  we  have  expressed  the 
radiation  integral  in  (2.12)  in  terms  of  a series  given  in  (2.24).  We 
have  shown  that  the  first  term  of  the  series,  viz.,  E^®^(u,v),  gives  the 
radiated  field  for  (u,v)  = (ug,VQ)  where  Ug,  Vg  is  typically  taken  to  be 
the  beam  maximum.  The  first  term  also  predicts  the  pattern  with  good 
accuracy  for  e not  too  close  to  Og  if  we  choose  s^  = 1,  since  the 
radiation  away  from  the  caustics  comes  mainly  from  the  rim  of  the  reflector 
(s  = 1)  and  the  particular  choice  of  s,  makes  a small  in  this  region  of  the 

a 

"effective  aperture."  For  the  intermediate  range  of  observation  angles, 
we  can  efficiently  compute  the  higher-order  integrals  in  (2.24),  i.e., 

Ep^,  etc.,  by  making  advantageous  use  of  the  recursion  relations 

of  the  Jacobi  polynomials.  In  fact,  among  other  reasons,  this  property 
of  the  Jacobi  polynomials  makes  them  uniquely  suited  for  our  purpose  of 
series  expansion. 
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In  the  next  section  we  present  some  numerical  results  illustrating 
the  application  of  the  theory  developed  in  this  paper. 

2.4  Steps  for  Numerical  Computations  and  Illustrative  Results 

Given  the  primary  pattern  of  a feed  and  the  parameters  of  the 
reflector  antenna,  the  numerical  computation  of  the  secondary  pattern  is 
carried  out  using  the  following  steps: 

(1)  The  current  distribution  3 on  the  reflector  surface  is  computed 
from  the  knowledge  of  the  feed  location  I and  its  primary  pattern  using 
the  physical  optics  approximation. 

(2)  f(s,i<i')  is  obtained  from  (2.4).  Note  that  p sin  e'  = r and 
r = as.  The  same  procedure  is  repeated  for  other  scalar  components  of  3 
(for  details  see  Appendix  II). 

(3)  Beam  deviation  factors  are  computed  using  the  method  in 
Appendix  I and  u^^,  v^^  are  derived  from  the  coordinates  of  the  phase 
center  of  the  displaced  feed.  We  also  calculate  6|^,  t|,,|  and  v^^  appearing 
in  (2.10). 

(4)  fj  is  derived  using  (2.11b). 

(5)  Jacobi  polynomials  are  generated  using  their  recursion 
relationship  [6]. 

(6)  The  two-dimensional  integration  represented  in  (2.21)  is 
performed  to  compute  C|j^  and  D^.  The  integration  may  be  carried  out  using 
any  of  the  standard  quadrature  routines,  e.g.,  the  Gauss  quadrature.  If 
a large  number  of  <))-coefficients,  i.e.,  n,  are  required,  one  may  perform 
the  (^-integration  via  the  FFT. 

(7)  The  coefficients  c"  and  are  inserted  into  the  recursion 

mm 

relation  (2.29)  to  generate  the  higher-order  coefficients  ^C„, 

(0)  (1)  (21 

(8)  Ej  , Ej  s . etc.,  are  successively  computed  using  these 
coefficients  and  Equations  (2.19),  (2.25)  and  (2.26)  are  added  up  to 
generate  Ej  and  subsequently  F. 

(9)  The  procedure  is  repeated  for  the  other  components  of  3 using 

the  corresponding  f.  | 

(10)  The  vector  far-field  pattern  is  constructed  by  appropriately  I 

summing  up  the  various  contributions.  | 


We  now  present  a few  numerical  results  to  illustrate  the  convergence  fj 

of  the  method.  Figure  2 shows  the  effect  of  increasing  the  number  of  | 


terms  in  the  Fourier  series  expansion  of  f,,  i.e.,  m , (=  M)  and  n 

1 IDd  X nid  X 

(=  N)  for  increasing  lateral  displacement  of  the  feed,  e.g.,  (or 
equivalently  kau^).  The  convergence  is  seen  to  be  quite  rapid  even  for 
several  beamwidths  of  scan  which  is  accompanied  by  a distortion  of  the 
secondary  pattern. 

Figure  3 shows  the  convergence  of  the  p-series,  i.e.,  (2.24).  It 
is  evident  that  even  for  a badly  distorted  pattern  only  one  term  of  the 
p-series,  i.e.,  is  sufficiently  accurate  for  computing  both  the 

amplitude  and  phase,  except  in  a relatively  small,  intermediate  region 
20  < kau  < 30.  However,  only  three  terms  give  sufficiently  accurate 
solutions  over  the  entire  range  of  observation  angles. 

Before  closing  this  section  we  remark  that  unlike  the  conventional 
methods  the  present  approach  is  very  well -suited  for  computing  radiation 
patterns  of  large  reflectors,  even  several  hundred  wavelengths  in  size 
and,  in  fact,  the  relative  advantages  of  the  series  methods  become  even 
more  evident  when  dealing  with  the  large  reflector  antennas.  Since  the 
p-series  becomes  more  rapidly  convergent  as  D/x  is  increased,  one  finds 


CONVERGED 

M-N-2 

- N - 3 

M- N -4 

M*N*S 


* ■* 


Fig.  2.  Series  convergence  as  a 
function  of  m and  n. 
f/D  = 0.5,  D/A  = 50. 


FIRST  TERM  OF  CXPONCNTiAL  SERIES  (f  - 1) 

2 TERMS  (p  • 2) 

3 TERMS,  (p- 2) CONVERGENCE  (EXACTI, 


Fig.  3.  Convergence  of  exponen 
tial  series  (2.24). 


that  the  usual  numerical  problems  of  storage,  loss  of  accuracy,  and 
astronomical  increase  of  computation  time  are  virtually  nonexistent  here. 
Furthermore,  Ej^^(u,v)  can  be  expressed  in  a scaled  form  such  that  the 
results  for  the  pattern  are  universal  and  independent  of  x.  These  added 
features  are  quite  useful  for  computation  and  display  of  the  radiation 
patterns. 


3.0  The  Offset  Reflector 


3. 1 Introduction 

The  offset  reflector  is  finding  increasing  use  in  many  communication 
satellite  antennas  because  of  its  freedom  from  blockage  effects  which  can 
affect  the  performance  of  a reflector  antenna  in  a deleterious  manner. 

The  analysis  of  offset  reflectors  is  made  difficult  by  the  lack  of  any 
possible  simplifications  due  to  symmetries  in  the  structure.  In  addition, 
cross-polarization  effects  and  pattern  distortion  due  to  scanning  become 
more  important  here  than  in  the  symmetric  case.  Thus,  accurate  but 
efficient  computation  of  the  vector  secondary  pattern  of  offset  reflec- 
tors is  highly  desired  in  many  applications. 

In  this  work  we  show  that  with  only  a relatively  small  additional 
effort  one  can  extend  the  series  method  of  Section  2,  developed  for  the 
symmetric  parabola,  to  the  offset  case.  We  proceed  to  demonstrate  this 
in  the  following  sections. 

3.2  Radiation  Integral  for  the  Offset  Parabolic  Reflector 

The  geometry  of  the  offset  reflector  together  with  the  location  of 
the  feed  plane  is  shown  in  Fig.  4.  Once  again  we  can  write  the  normalized 
radiation  integral  in  terms  of  the  induced  current  j on  the  surface  of  the 
reflector  as 

, . . r r -jk(p'-p-L) 

£(0,4.)  = (!  - RqRq)  J J j e ^ ds  (3.1) 

SURFACE 

which  is  identical  to  (2.1).  Recall  that  the  optical  phase  of  j has  been 
factored  out  explicitly  relative  to  an  ideal  phase  center  located  at  e 
(see  Fig.  4).  Also,  the  integral  is  defined  on  the  surface  A of  the 
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Fig.  4.  Geometry  of  offset  parabola 


reflector  and  not  in  terms  of  an  approximate,  ray  optically  derived  field 
in  the  projected  aperture. 


Introducing  a mathematical  transformation,  we  can  rewrite  the  integral 
in  terms  of  coordinates  r,  ({i'  of  the  projected  aperture  of  the  parent 
parabola  (see  Fig.  4).  We  get,  for  a typical  scalar  component  of  E,  say, 
F(0,(t>),  integrals  of  the  type 


where  f is  related  to  the  corresponding  scalar  component  of  J via  the 

Jacobian  J of  the  transformation  between  the  parabolic  surface  and  the 

1/2 

projected  aperture.  The  Jacobian  J is  given  by  (1  + r/2f) 


Using  the  small  displacement  approximation,  we  can  write 


■cr^~ 


As  mentioned  earlier,  this  approximation  is  easily  generalized,  if  desired, 
and  is  not  really  necessary  for  our  work.  However,  the  small  displacement 
approximation  is  often  used  because  of  its  simplicity  and  because  it  is 
sufficiently  accurate  as  long  as  the  feed  displacement  is  not  too  large. 

For  a parabola, 

p = 2f  + p cos  e'  and  sin  e'  = r/fM  (3.4) 

with  M = 1 + r^/4f^  , 

and,  hence, 

p*  = 2f  + p cos  0'  - p • e 

= 2f  + p cos  0‘  - [(e  r/fM)  cos  (}>'  + {e„r/fM)  sin  (|)']  . (3.5) 

A y 

Recognizing  that 

p-RQ  = psine'sine  cos  (<|)  - <!>')+  p cos  e cos  e' 

= r sin  0 cos  {(()-(}>')+  p cos  0 cos  0’  , (3.6) 

we  can  rewrite  (3.1)  as 


F(e,^)  e^^^f  ^ [f  g-jkpcos0'(l-cose)  ^^^^z 


jke,COS0 ' 


where 


and 


X e 


-jkSjr(u^cost'*v^sfn,') 

X 0 


r dr  d((>' 


= 


- (r/2f)r 


1 + (r/2f)" 

.1  (v  + V ) r p ^ 

“ " Tu  +"u^)  i = /(v  + + (u  + u^)^ 

G e 


(3.7) 

(3.8) 

(3.9) 
(3.10) 


Examining  (3.7),  we  note  the  following  features  of  the  integrand. 

In  addition  to  the  function  f,  which  is  realted  to  the  induced  surface 
current  on  the  reflector,  the  integral  has  four  exponentials.  The  first 
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of  these  is  a function  of  both  the  observation  angle  e and  integration 
variable  e‘(r).  The  second  and  third  exponents  depend  on  the  integration 
variables  r and  <{>'  only.  In  the  second  exponent  we  have  deliberately 
left  to  be  nonzero,  since,  in  contrast  to  the  symmetric  reflector, 
the  plane  of  displacement  of  the  feed  for  minimum  distortion  is  not 
entirely  transverse  for  the  offset  case.  We  will  soon  demonstrate  this 
point  when  we  derive  a formula  for  the  optimum  feed  displacement  plane. 

The  third  exponent  contains  u^  and  v^,  which  define  the  lateral  coordinates 

of  the  displaced  feed.  Finally,  the  last  exponent  is  the  familiar  Fourier 
transform  kernel  expressed  in  cylindrical  coordinates. 

It  is  desirable  for  later  manipulation  to  shift  the  origin  of  the 
integration  coordinates  to  the  center  of  the  projected  aperture.  To  this 
end  we  introduce  the  transformation  of  coordinates  from  (r,(()')  to  (rj^,(()p 
system  (see  Fig.  4)  via  the  equations 

r cos  <()'  = d^  + r^  cos  41^  (3.11a) 

r sin  4.'  = r^  sin  41^  (3.11b) 

Writing  F(u,v)  = E^(u,v},  we  get  from  (3.7), 
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0 JO 
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X e e 

jkrir2Cos(a-4)2) 


X e 
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(3.12) 


with 


^(^2,42)  = f(r,4') 
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(3.13) 


Note  also  that  the  r2  and  r are  related  by 


r^  = r2  + d2  + 2r2d2  cos  4*2 


(3.14) 
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At  this  point  it  is  worthwhile  to  point  out  some  of  the  similarities  and 
differences  between  the  radiation  integrals  for  the  symmetric  and  offset 
reflectors.  For  both  of  these  cases  the  radiation  integral  reduces  to  the 
Fourier  transform  form  if  the  first  exponent  is  small,  which  occurs  when 
e = 0.  In  the  symmetric  case  the  function  z(rj,ij)^)  is  independent  of  the 
azimuthal  coordinate  and,  in  addition,  is  a quadratic  in  the  radial 
coordinate.  This  special  property  of  the  z-coordinate  of  a symmetric 
parabola  allowed  us  to  express  the  radiation  integral  in  the  form  of  an 
analytically  continuous  series  whose  higher-order  coefficients  are 
derivable  from  those  of  the  corresponding  zero-order  term  via  the  re- 
cursion relationship  of  the  Jacobi  polynomials  given  in  the  previous 
section. 

It  would  appear  at  first  sight  that  this  unique  and  advantageous 
feature  is  no  longer  available  for  the  offset  case.  Since  the  function 
z(r^,(jij^)  is  neither  independent  of  the  azimuthal  variable  nor 
is  it  a quadratic  function  of  r^.  However,  we  show  in  Sec.  3.4  that 
by  using  a suitable  transformation  of  the  observation  variables  n and  a 
we  can  still  recast  (3.12)  into  a form  which  has  all  of  the  desired 
characteristics  of  the  symmetric  case  including  the  numerical  efficiency 
and  ease  of  computation. 

Before  discussing  this  topic,  however,  we  will  briefly  digress  to 
consider  the  problem  of  choosing  the  optimum  scan  plane,  i.e.,  as  a 

function  of  (e  ,e  ),  for  offset  reflectors. 

X y 

3. 3 Choice  of  Optimum  Scan  Plane 

We  return  to  (3.12)  and  rewrite  the  third  exponent,  which  contains 
as 

jk(6^d^Ug  - cos  0')  = jkej[djU^  - 2c^]  + jkc^  (3.15) 

where  was  defined  in  (3.8).  This  allows  us  to  rewrite  (3.12)  as 


-jk(djUj-e^)-jkdjU 


Ej(u,v)  = Ej(u,v)  e 

r^l  r2v  r.  -jkz(r,,((>,)[l-cos0] 

= L L ' 


-jk6jr^tgC0s(Y  -<>j)  -jkej[d^Ug-2ej,] 

X e e 

jknr,cos(a-(|),)'> 

X e ^ ^ J (3.16) 

where  we  have  used  u = n cos  a.  It  may  be  worthwhile  to  point  out  that 
-jkd.u 

e ^ represents  the  phase  factor  introduced  in  the  expression  for  the 
radiation  integral  due  to  a shift  of  the  origin  from  the  focus  of  the 
parent  parabola  to  the  center  of  the  projected  aperture. 

The  criterion  we  choose  for  determining  EylEy.Cw)  is  that  an  optimum 

z X y 

choice  should  minimize  the  phase  error  terms  represented  by  the  second 
and  third  exponential  factors  in  (3.16),  exclusive  of  linear  phase  tapers 
which  may  contribute  to  beam  shift.  Explicitly,  we  deal  with  the 
exponent  h(rj,<j)j) 


h(r^,(t)p  = k6^[djU^  - 2e^]  + kBj^lr^t^  cos  (y^  - 4>^)] 


(3.17) 


We  note  that  the  first  term  in  the  square  bracket,  i.e.,  (djU^  - 2e^] , is 
independent  of  r^  and  (ji^,  whereas  the  second  term  contains  the  factor 
r^^tg  cos  (yj  - i))j),  which  has  a single  cosinusoidal  variation  as  a function 
of  the  azimuthal  angle  Since  both  of  these  have  a common  multiplicative 
factor  kBp  one  might  be  tempted  to  conclude  that  the  minimum  distortion 
due  to  h(rp(j)j)  results  when  we  set 

= djU^  . (3.18) 

However,  as  shown  below,  the  correct  choice  differs  from  that  given  in 
(3.18)  because  the  behavior  of  Bj  as  a function  of  (r^ifi^)  can  not  be 
ignored.  To  this  end  we  examine  B^  and  express  it  as 
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Si  = 


(r/ZfV 


(r^/2f)2  + (dj/2f)^  + (rj/2f)(dj/2f)  cos 


^ 1 + (r/2f)^  1 + (rj/2f)^  + (dj/2f)^  + (rj/2f )(dj/2f)  cos 


(3.19) 


which  may  be  approximated  to  give 


Q(r^  + dj)  + 2rjdj  cos  (|)0 


(3.20) 


We  therefore  have 


h(r^,(})P  = 


+ r 


+ r' 


djtg  cos  (Yj  - <tij)  + 2d^(d^u^  - 2e^)  cos  ())^ 

(djtg  cos  Yg  + djUg  - 2e^)  + djtg  cos  (y^  - 2(|.^)j 


+ r^t^  cos  (Yg  - <^j) 


(3.21) 


We  note  in  (3.21)  that  the  parameter  appears  only  at  three  places. 

The  first  term  containing  is  a constant  in  r^  and  merely  introduces  a 
constant  phase  factor  in  the  radiation  integral  without  contributing  to 
distortion.  The  second  term,  linear  in  r^,  also  contains  e^.  However, 
this  exponential  factor  only  contains  cos  (jij^  and  sin  ()ij  types  of  terms, 
and  hence  can  be  straightforwardly  combined  with  the  Fourier  transform 
kernel  -exp[jkn  cos  (y^  - to  derive  the  direction  of  the  shifted  beam. 
We  show  this  explicitly  in  the  next  section.  Meanwhile,  we  remark  here 
that  because  of  its  special  form  the  linear  term  in  r^  does  not  introduce 
any  distortion  terms  either;  therefore,  we  turn  our  attention  to  the 
quadratic  term  in  rj^.  It  is  evident  that  the  first  term  inside  the  bracket 
is  independent  of  (|ij,  whereas  the  second  term  has  a cyclical  variation 
along  the  azimuthal  direction.  Consequently,  the  distortion  in  h(rj,(|>j) 
as  a function  of  is  minimized  by  enforcing  the  condition 

*^l^s  ^s  ^ *^l^s  ” ^^z  ~ ^ (3.22a) 
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or 

= (l/2)[djtg  cos  Ys  + djUg]  (3.22b) 

which  can  be  further  simplified  to  give 


e 


z 


(3.23) 


by  recognizing  that  t^  cos  = u^.  Comparison  of  (3.23)  with  (3.18)  shows 
that  the  two  expressions  differ  by  a factor  of  two.  We  now  show  that 
(3.23)  implies  that  the  optimum  plane  of  displacement  of  the  feed  is 
approximately  the  plane  normal  to  the  line  OP  where  P is  the  point  on  the 
surface  of  the  reflector  whose  projection  in  the  focal  plane  coincides  with 
the  origin  (x  = d^,  r^  = 0)  of  the  projected  aperture. 

Along  this  plane  and  satisfy  (see  Fig.  5) 
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Fig.  5.  Geometry  of  optimum  plane 
for  feed  displacement. 


for  f/D  in  the  range  0.25  < (f/D)  < 1, 
since  typically  (dj/2f)^  < {l/[8(f/D)]} 
for  an  offset  reflector. 

The  result  shown  in  (3.25),  though 
simple,  is  quite  important  from  a 
practical  point  of  view.  In  the  past, 
the  choice  of  the  optimum  scan  plane 
has  been  based  on  experimental  trials 
and  a systematic  derivation  of  the 
equation  for  the  plane  has  not  been 
available.  Before  closing  this  sec- 
tion we  add  the  remark  that  for  the 
symmetric  case,  dj  = 0;  hence,  the  op- 
timum scan  plane  reduces  to  the  focal 
plane  e = 0. 


25 


3.4  Series  Development  of  the  Radiation  Integral 

We  now  return  to  (3.16)  and  examine  the  possibility  of  developing  a 
series  expansion  for  this  radiation  integral  in  the  same  manner  as  was 
done  for  the  symmetric  case.  We  note  first  of  all  that  in  the  absence 
of  the  first  exponential  factor  exp[-jk2(l  - cos  e)]  the  integral  (3.16) 
would  reduce  to  a Fourier  transform  form  and  could  therefore  be  efficiently 
evaluated  by  any  number  of  available  methods.  Next,  we  recall  that  this 
factor  was  also  present  in  the  symmetric  case  and  that  we  were  able  to 
handle  it  by  using  the  recursion  relations  of  the  Jacobi  polynomials 
because  the  z(r,(ti)  for  the  symmetric  reflector  assumed  a very  special  form, 
viz.,  quadratic  in  r and  independent  of  (p.  The  apparent  difficulty  in 
using  the  same  procedure  for  the  offset  case  stems  from  the  fact  that  z is 
now  a function  of  both  the  integration  variables  r^  and  4>j,  and  the 
recursion  relationships  of  the  Jacobi  polynomials  cannot  be  directly  used 
to  efficiently  evaluate  the  coefficients  of  the  Jacobi  polynomial  series 
for  the  higher-order  terms  of  the  exponential  series  derived  from  the 
expansion  of  exp[-jkz(l  - cos  e)].  It  is  extremely  fortuitous,  however, 
that  certain  manipulations  of  the  integrand  in  (3.16)  allow  us  to  retain 
all  of  the  advantages  of  the  Jacobi  polynomial  series  without  any  penalty 
whatsoever.  We  proceed  now  to  elaborate  on  this  point  in  the  following 
paragraphs. 

As  a first  step  we  explicitly  write 


z(rj,it>j)  = p cos  9' 


2f  cos  e' 

1 - cos  e' 

2 2 

r^  + dj  + 2rjdj  cos 
4f 


4r 


(3.26) 


for  the  offset  parabola.  Thus 

d^  - 4f^  r.d 

z(rj,(j>j)(l  - cos  0)  = (1  - cos  0)  + cos  (1  - cos  0) 

+ ^ (1  - COS  9)  . (3.27) 


We  now  note  the  following  features  of  the  three  terms  appearing  in  the 
r.h.s.  of  (3.27).  The  first  term  is  independent  of  the  integration 


c- 


variables  and  depends  only  on  the  observation  angle  e.  The  second  term 
has  a special  form  in  that  it  is  linear  in  r^  and  has  a cos  variation. 
Finally,  the  last  term  is  quadratic  in  r^  and  independent  of  We  will 
soon  make  advantageous  use  of  these  properties  of  the  first  exponent  in 
(3.16). 

Turning  now  to  the  second  exponent  in  (3.16),  we  write 

Bi^its  cos  (yj  - ((>j)  = 32rj(Ug  cos  ((.j  + v^  sin  (f.^) 

= + C^j)  cos  ({>2  + (Bj^s  + Cy)  sin 

- rjlC^  cos  ({.j  + sin  (3.28) 

where  C^,  are  constants  yet  to  be  specified. 

Finally,  we  rewrite  the  term  nr^  cos  (a  - ^i^)  in  the  last  exponent  as 

nrj  cos  (a  - (tij)  = rjln  cos  a cos  (l)^  + n sin  a sin  (ji^] 

= r2[(u  + Ug)  cos  + (v  + v^)  sin  (ji^]  . (3.29) 

Using  (3.27),  (3.28)  and  (3.29),  we  can  now  combine  certain  terms  from  the 
various  exponents  in  (3.16)  that  contain  the  factors  r^  cos  and 
r^  sin  (i)^,  respectively.  We  choose  to  collect 

r^  cos  (jij  terms 
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r^  sin  (ji^  terms 
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from  the  various  exponents  and  introduce  a new  pair  of  observation  vari- 
ables n2  3nd  02  via  the  definitions 


02  cos  “1  = “ ^ (1  - cos  0)  + + (u  + Ug) 


^2  sin  a2  = Cy  + (v  + v^) 


(3.29a) 

(3.29b) 
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hence, 

rj[nj  cos  aj  cos  (|ij  + sin  sin  41^]  = cos  (a^  - 41^)  . (3,30) 


Using  (3.27)  through  (3.30)  in  (3.16),  we  can  rewrite 
jk.Qdj-4f2)/450-cose)  j(La"'^ao) 
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(3.32) 

(3.33) 

(3.34a) 

(3.34b) 

(3.o4c) 

(3.34d) 


Comparison  of  (3.31)  with  (2.11a),  (3.32)  with  (2.11b),  (3.33)  with 
(2.9d)  and  (3.34)  with  (2.7)  reveals  an  almost  one-to-one  correspondence 
between  the  equations  for  the  radiation  integral  for  these  two  cases, 
provided  that  the  observation  coordinates  (n.a)  for  the  symmetric  case  are 
replaced  by  (n^.a^  for  the  offset  case  and  certain  suitable  proportion- 
ality factors  are  introduced  in  defining  and  ^2- 

It  is  evident  that  the  processing  of  (3.31)  can  follow  exactly  the 
same  treatment  as  (2.11).  That  is,  the  function  fj(s,4>j)  can  again  be 
expressed  in  tprms  of  a two-dimensional  Fourier  series  in  s and  4>j  in  the 
same  manner  as  in  (2.13)  and  the  exponential  series  for  (3.31)  be  obtained 
by  expanding  exp[-jA]  as  before.  The  evaluation  of  E2^^(u,v),  as  well  as 
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of  the  higher-order  terms  of  the  series,  viz.,  etc.,  defined  in 

the  same  manner  as  in  (2.19),  (2.24),  (2.25)  and  (2.26),  also  proceeds  in 
an  identical  manner  with  the  use  of  the  recursion  relation  for  the  Jacobi 
polynomials  for  the  higher-order  terms.  In  short,  all  of  the  advantages 
of  the  Jacobi  polynomial  method  of  evaluating  the  radiation  integral  are 
maintained  in  the  offset  case  without  any  penalty. 

Before  closing  this  section,  it  may  be  worthwhile  to  make  some 
comments  on  the  new  observation  variables  n^,  oj  defined  in  (3.29).  We 
note  first  of  all  that  the  origin  of  the  observation  coordinates,  say, 

Uq,  Vq,  may  be  fixed  by  letting  = (dj/2f)(l  - cos  Og)  - (ug  + u^)  and 
Cy  = -(v  + v^),  since  nj(ug,Vg)  then  becomes  identically  zero.  Typically, 
though  not  necessarily,  (ug,Vg)  are  chosen  to  coincide  with  the  beam 
maximum,  in  which  case 

Un  = , Vo  = v„ 

0 m 0 m 

Cnn,  “ (d./2f)(l  - cos  0„)  - (u„  + u.) 
urn  1 m'  ' m s 

vm  ' m s' 
and  we  have  from  (3.29) 

d 

cos  aj  = -^  (cos  6 - cos  0^)  + (u  - u^^) 
ni  sin  aj  = V - v^ 

which  are  somewhat  more  convenient  to  use  than  (3.29a)  and  (3.29b). 

In  order  to  determine  u^  and  v^  for  a given  u^  (=  (= 

it  is  necessary  to  know  the  beam  deviation  factors  and  B^.  The 
beam  deviation  factors  for  an  offset  parabola  can  again  be  expressed 
in  terms  of  the  primary  feed  pattern  [7]. 

3.5  Steps  for  Numerical  Computation  and  Illustrative  Results 

The  steps  for  numerical  computation  of  the  secondary  pattern  for  a 
given  primary  pattern  of  the  feed  follow  pretty  much  along  the  lines 
described  in  Sec.  2.4. 

Some  illustrative  results  of  offset  pattern  computation  are  shown 
in  Fig.  6.  Rapid  convergence  is  achieved  (see  Fig.  6b)  even  when  the  beam 
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is  considerably  distorted  as  compared  to  the  focal-feed  pattern  (Fig.  6a), 
which  of  course  converges  even  faster.  Once  again  we  mention  that  the 
computation  time  for  secondary  patterns  is  much  more  rapid  here  as  compared 
to  the  conventional  numerical  integration  or  the  spherical  wave  expansion 
and  that  the  relative  advantage  becomes  even  more  apparent  as  the  dish 
size  gets  larger. 

4.0  Contour  Beam  Synthesis  Using  Reflector  Antennas 
4. 1 Introduction 

Reflector  antennas  typically  radiate  narrow  pencil  beams  and  even 
though  some  control  over  the  shape  of  the  radiation  pattern  can  be 
exercised  by  perturbing  the  reflector  or  subreflector  shapes  and  their 
contours,  this  step  alone  seldom  offers  the  flexibility  needed  to  closely 
match  the  coverage  to  a given  geographical  area,  e.g.,  the  Eastern  time 
zone  (ETZ)  of  the  U.S. 

A much  more  satisfactory  approach  to  synthesizing  contour  beams 
entails  a superposition  of  a multitude  of  pencil  beams  generated  by  a 
cluster  of  feeds  whose  complex  excitation  coefficients  are  adjusted  to 
obtain  the  desired  secondary  power  pattern.  In  this  section  we  address 
ourselves  to  this  problem  of  designing  cluster  feeds  for  reflector 
antennas  with  the  objective  of  synthesizing  contour  beams.* 

Conventionally,  the  synthesis  problem  is  attacked  via  the  so-called 
spotlighting  technique  in  which  the  complex  amplitude  of  the  feed  that 
generates  a beam  falling  within  the  contour  is  set  equal  to  1,  and 
identically  zero  otherwise.  Although,  generally  speaking,  this  approach 
is  often  quite  satisfactory,  it  nevertheless  suffers  from  the  drawback  that 
there  is  very  little  control  over  the  slope  of  the  contour,  sidelobe  level, 
etc.  Perhaps  even  more  important  is  the  fact  that  the  method  does  not 
take  into  account  the  effect  of  beam  distortion  with  lateral  displacement 
of  the  feed,  which  may  be  substantial  even  when  the  beam  scan  is  only  four 
or  five  beamwidths. 

Neither  does  the  spotlighting  method  necessarily  generate  an  optimum 
or  "best-fit"  solution  to  the  specified  contour  pattern,  which  is  typically 
a power  pattern,  since  it  attempts  to  approximate,  at  best,  the  specified 
field  pattern  while  arbitrarily  assigning  the  phase  distribution  within  the 


* The  material  is  based  on  a forthcoming  paner  [8]. 


contour  to  be  uniform.  Fortunately,  however,  the  spotlighting  techniques 
can  still  be  used  to  generate  results  that  are  very  close  to  the  "optimum" 
one  — typically  in  an  extremely  efficient  manner,  providing  the  constituent 
beams  remain  relatively  undistorted  as  these  beams  are  scanned.  Obviously, 
this  is  not  possible  to  achieve  with  a single  feed  unless  the  reflector  is 
shaped  to  reduce  the  scan  distortion,  a task  that  is  not  easily 
accomplished,  and  in  fact,  no  systematic  means  for  doing  it  is  available 
at  this  time. 

Nevertheless,  as  we  show  in  Section  4.2,  a cluster  of  seven  feeds  can 
be  synthesized  for  each  scan  position  such  that  the  beam  distortion  with 
lateral  displacement  is  substantially  reduced.  Superposition  of  these 
composite  beams  using  the  sootlighting  technique  can  then  be  employed  to 
synthesize  a desired  contour  beam. 

As  mentioned  earlier,  the  optimization  approach  has  both  the  flexibil- 
ity and  capability  to  handle  the  general  synthesis  problem  since  it  can 
deal  directly  not  only  with  the  specified  power  pattern,  but  can  also  be 
programmed  to  achieve  some  control  over  the  ripple,  slope  and  sidelobes 
associated  with  the  contour.  We  discuss  this  method  in  Section  4.3  and 
describe  certain  algorithms  for  extracting  a numerical  solution  to  the 
optimization  problem.  Finally,  in  Section  4.4,  we  present  some  examples 
of  application  of  both  of  these  methods  to  the  problem  of  synthesizing  the 
ETZ  contour. 

4.2  Synthesis  of  Cluster  Feeds  to  Reduce  Pattern  Distortion 

In  this  section  we  discuss  two  approaches  for  reducing  the  distortion 
in  the  secondary  pattern  due  to  laterally  displaced  feeds.  We  begin  with 
the  dominant  term  for  the  radiation  integral,  viz.,  E^^^(u,v),  which  is 
given  by  [see  (2. lial] 

(•1  -Z-Ft 

En(‘J«v)  = ! s ds  ' 6^'  f.{s,4>')  exp[^  jkDsn  cos  ((|>'  - a)]  . (4.1) 

^ -1 0 -0  ^ 

We  restrict  ourselves  to  the  case  of  symmetric  reflectors  and  remind  the 
reader  that  the  derivation  of  (4.1)  was  given  in  Section  2.  For  convenience, 
we  repeat  some  of  the  pertinent  definitions: 


= r'/a  = Zr'/D 


(4.2a) 


Eu  - + (v  - 

tan"^  Q(v  - Vj^)/(u  - Uj^^)3 


(4.2b) 


(4.2c) 


= sin  0 cos  <ti  , V = sin  0 sin 


(4. 2d) 


The  location  of  the  beam  maximum  is  given  by  Uj^^,  v^^  defined  as 


u = U|v,  = B^j(e/f)  cos  Y. 


(4.3a) 


V = = By(e/f)  sin  Yg 


(4.3b) 


where  B^  are  the  beam  deviation  factors.  The  function  is  derived 


from  the  primary  feed  pattern,  and  the  location  of  the  feed  was  defined 
in  (2.11b).  It  can  be  rewritten  as 


fj(s,V)  = P*(s)  Q*(s.<^')  f(s) 


(4.4) 


where  f(s)  is  associated  with  the  illumination  function  of  a focal  feed 
[see  (2.4)]  and  is  assumed  here  to  be  circularly  symmetric.  The  functions 
P(s)  and  Q(s)  are  given  by 


P*(s)  = exp 


2 

-2  ke(e/f)  ^ ■%  - kf(l  - N)  1 - 

1 + N L 


1 - (B^e/f) 


Q*(s,4>')  = exp{-jt0e  COS  (<[.'  - Yf..)]} 


(4.5a) 

(4.5b) 


where 


^ " 4(w] 

V,-  — ^ 


(4.6a) 


6 = s-i  ks(f/D)“^  B„  + (1  + 


(4.6b) 


Y,,  = Y If 

'M  's 


(4.6c) 


Note  that  P and  Q are  distortion  factors  introduced  by  the  displacement  of 


the  feed,  i.e. , due  to  nonzero  e and  e . 

X y 


The  secondary  pattern,  which  can  be  calculated  using  the  series  expan- 
sion method  described  in  Section  2.1,  is  found  to  have  a distortion  for  wide 
angle  scan.  We  have  calculated  the  secondary  pattern  with  the  following 
choice  of  parameters 


D = 50  A , f/D  = 0.5 


(4.7) 


f(s)  = (-cos  6 


(4.8) 


where  M was  defined  in  (4.6a)  and  cos  0'  is  related  to  s by  the  equation 
appearing  underneath  (2.10).  The  exponent  n controls  the  edge  taper  of 
the  illumination,  e.g.,  n = 2.2538  for  a -10  dB  taper. 


The  secondary  pattern  in  dB  is  plotted  in  Fig.  7 as  a function  of  kan. 
We  observe  that  the  beamwidth  is  approximately  2 x 1.8  in  the  ka  space  and 
that  the  sidelobes  decrease  monotonically  from  the  -26  dB  level  of  the 


first  one. 


In  Fig.  8 we  have  shown  the 
results  for  the  secondary  pattern 


obtained  with  the  feed  displaced  from 
the  focus  to  e = 3a,  = tan"^ 


= IT.  We  observe  from  Fig.  8 that  the 
beam  becomes  quite  broad  in  the  scan 
plane  due  to  phase  errors  in  the 
aperture.  This  error  also  creates 
higher  sidelobes,  e.g.,  -13  dB  for 
the  case  shown.  The  effect  on  the 
pattern  in  the  plane  orthogonal  to 
the  displacement  of  the  feed  is 
rather  minimal . 


We  now  address  ourselves  to  the 


Secondary  pattern  of  a 
parabolic  reflector  due 
to  a focal  feed. 

D/X  = 50,  f/D  = 0.5 
(10  dB  edge  taper). 


problem  of  compensating  for  the  phase 
error  in  the  aperture  introduced  by 
the  displacement  of  the  feed.  For 
this  purpose,  we  consider  an  array 
of  seven  feed  elements  distributed 
in  a hexagonal  lattice  shown  in 
Fig.  9.  Our  problem  is  to  choose 
the  excitation  coefficients  for  the 


seven  feed  elements  such  that  the 


SINGLE 

ELEMENT 


Fig.  8.  Secondary  pattern  distortion  due  to  a displaced  feed.  Parameters 
same  as  in  Fig.  7.  Note  on  focus  pattern  shown  in  dotted  line. 

distortion  factors  are  compensated 
for  as  much  as  possible.  The  function 
for  the  cluster  of  seven  elements 
takes  the  form 

f^(s,(t,‘)  = P*(s)Q*(s,<fr')f(s)A(s,<(.') 

(4.9) 

where  A(s,(j)')  is  the  array  factor 
given  by 

A(s,<(>') 

= 1 + I I exp/j<d 
m=2  L 

and  where 

Fig.  9.  Feed  as  a periodic  array  of 

equilateral  triangular  < = k sin  9'  = ks/ 

spacing.  Elements  1 to  7 
form  a typical  cluster.  (4.10b) 


y 
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Ijj,i  ni  = 2, 3, 4,..., 7 in  (4a)  are  the  cotnple  excitation  coefficients  of  the 
ring  elements  and  without  loss  of  generality,  we  have  set  the  coefficient 
of  the  center  element  equal  to  1. 

It  is  evident  from  (4.9)  that  to  eliminate  the  distortion  from  f^  our 
goal  is  to  adjust  A(s,(t)')  such  that,  ideally, 

A(s,(t)')  = [P*Q*]“^  , P.Q  = 0 < s < 1,  0 < (<.’  < 2it  . (4.11) 

Strictly  speaking,  it  is  not  possible  to  achieve  tfie  condition  enunciated 
by  (4.11)  exactly  since  A has  only  twelve  degrees  of  freedom.  However, 
we  can  satisfy  (4.11)  approximately  by  matching  the  leading  Fourier 
coefficients  of  both  sides  of  this  equation.  We  discuss  a method  for 
doing  this  in  the  following  paragraphs. 

The  problem  of  expressing  the  Fourier  coefficients  for  the  expansion 
of  A(s,(j)')  is  simplified  considerably  if  instead  of  working  directly  with 
the  unknown  excitation  coefficients,  we  introduce  a new  set  of  unknowns 
via  the  equation 

= I S , m = 2,3,4,. ..,7  (4.12) 

q=0  ^ 

which  is  an  orthogonal  transformation  of  the  Substituting  (4.12)  into 
(4.10)  gives 


A(s,V)  = 1+1  S B (s,4.‘) 

q=0  ^ ^ 


(4.13a) 


where 


^(s,(|)')  = exp^’Kd  cosQ'  - (m  + 1)  1]]+  jq(m  + 1)  ^ • (4.13b) 


Making  use  of  the  identities 


^jxcose  . J ♦ e-j"] 


m=0 

5 „ (6,ifn  + q = 6fc 

I exp  j(q  + n)  = 

_r»  J 


n=0 


0,  otherwise 


(4.14a) 

(4.14b) 
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in  (4. 13),  we  get 


Bq(s.<fr')  = 6 I (4.15) 

and  eventually  the  Fourier  series  for  A(s,(|)') 

A(s,(^')  = 1 + 65q  I e^(j)®^  Jg^(Kd)  cos  6t4.' 

£ = -00 

♦ 6 i (J)“*2  } Wi)L  eJ(6‘+2)*'  , 5^ 

£=-00 


At  this  point,  we  note  that  the  product  P.Q,  which  we  are  attempting  to 
approximately  equate  to  A(s,()>'),  is  an  even  function  of  (<(i'  - Y|^^),  whereas 
the  Fourier  series  given  in  (4.16)  is  not.  In  fact,  it  turns  out  that  with 
the  limited  number  of  degrees  of  freedom  in  our  hands,  it  is  not  possible 
to  make  the  function  A(s,((>')  even  in  (a*  - i'or  an  arbitrary  Y|v^* 

However,  it  is  possible  to  enforce  the  condition  of  evenness  to  the  five 
Fourier  coefficients  of  for  p = 0,±1  and  ±2.  We  can  show  that  this 

is  possible  by  imposing  the  two  constraints 


(4.17) 


which  reduces  the  number  of  unknown  coefficients  from  six  to  four.  Denoting 
these  new  unknowns  S . we  relate  them  to  S via 


We  can  also  relate  the  to  the  original  unknowns  as  follows 


In,  = S 
m 


^0  ^ Sq  cos  q + 1)  I - 

+ $2  cos  sQiti  + 1)  j , m = 2, 3, .,.,7 


(4.19) 


The  Fourier  series  for  the  array  factor  becomes 


A(s.^‘)  = 1 + 6Sq  I cos 


+ i2s^  I + 1)  t'  - Y^] 

l=-CO 

+ 12S2  I + 2)  - 2y„1 


(4.20) 


+ 6S3  I (j)®^^^  [(6£  + 3)  <1.*  - 3y,^]  . 

Jl  = -00 

We  note  from  (4.20)  that  only  the  cos  [p((|)'  - Yj^fj)]  terms  appear  in  the 
series  for  p = 0,  1 and  2,  thus  ensuring  the  symmetry  we  were  seeking  in 
A(s,(j)')  up  to  these  harmonics. 

Next  we  turn  to  the  problem  of  determining  by  performing  an 
approximate  matching  of  the  array  factor  A to  the  distortion  factor  P . Q. 
We  show  that  by  matching  the  appropriate  Fourier  coefficients  of  A and 
P . 0 we  can  derive  closed  form  solutions  for  the  coefficients  S^. 

Let  us  first  consider  the  case  of  Y|^  = 0.  The  function  P . Q may  be 
written  in  a Fourier  series 


P . Q = P(s)  z J Qn(s)  cos  n(j)' 


Q„(s)  . v.„(6.) 


(4.21) 


From  (4.20)  we  can  write 


A(s,((i' ) = E j A^(s)  cos  n({.' 


t 

\ 


■PPP^S 


with 


Aq(s)  = 1 + 6SnJn(<d)  , and 


0“0' 


A (s)  = 12S| 


n+6Ji.  n 


J.M) 


n = 1,2 


(4.22a) 

(4.22b) 


Equating  the  Fourier  coefficients  we  get 


A^(s)  = P(s)  Q^(s)  , for  0 < s < 1 and  n = 0,1,2,. 


. (4.23) 


It  is  evident,  once  again,  that  with  only  four  unknown  coefficients  it  is 
impossible  to  satisfy  (4.23)  for  all  s.  The  best  we  can  hope  for  is  a 
minimization  type  of  solution,  e.g.,  one  that  minimizes  the  error*  A defined 
by 


A = 


sds  W(s)  I lA  (s)  - P(s)  Q_(s)l‘ 
0 n=0  " " 


(4.24) 


In  general,  we  might  expect  (4.24)  to  lead  to  a 4 x 4 matrix  equati 


on 


to  be  solved  for  S^.  Fortunately,  however,  it  is  possible  to  diagonalize 


the  resulting  equation,  enabling  one  to  extract  the  solution  without  the 
need  for  matrix  inversion.  To  develop  this  diagonalized  form,  we  first 
write  A as 


q=0 


(4.25) 


where 


oo 

Aq  = (sds)  W(s)  |6°  + 6SQjg^(<d)  - P(s)  Jgn^ee)!^  (4.26a) 


1 


BESl  /, 


rcopY 


We  note  that  depends  only  on  one  S,  viz.,  S^.  Thus  the  individual  S^'s 
can  be  determined  by  minimizing  the  corresponding  A^'s,  one  at  a time.  By 
requiring 


— £ = 0 = 0 
3S'  ^ ’ 3S"  ^ 

q q 


(4.27) 


where  - jS|j,  we  can  derive  the  closed-form  solutions 


Sq- 


2 oo  ^ 

(sds)  W(s)  I 
0 n=0 


S = ^ 


p(s)  Jgn(ee) 


sds  W(s)  I [JfijKd)]^ 
0 n=0 


(sds)  W(s)  P(s)  I 

n = -a>  'w. 


'^6n+q^'^‘^^  '^en+q^^^^ 


sds  W(s)  I [Jgn+q^'"'^)^ 


0 n=-” 

q = 1,2.3  . 


(4.28a) 


(4.28b) 


Equation  (4.28)  represents  the  solution  for  we  were  seeking.  We  can  of 
course  determine  the  I^^'s  from  (4.19)  once  are  known. 

The  generalization  of  the  solution  to  an  arbitrary  y^/i  will  now  be 
discussed.  First,  if  Y|v]  = nu/G,  where  n is  a positive  or  negative  integer, 
the  entire  procedure  described  above  for  = 0 goes  right  through  without 
change.  For  Yf^i  f nir/e,  it  is  no  longer  possible  to  match  all  of  the 
(j)-harmonics  but  only  p = 0,  ±1  and  ±2  because  we  do  not  have  a sufficient 
number  of  degrees  of  freedom.  Once  this  fact  is  accepted,  the  resulting 
equations  and  the  solutions  for  remain  unchanged  from  (4.26). 

4. 3 Illustrative  Results  of  Cluster  and  Contour  Beam  Synthesis 

At  this  point  it  will  be  worthwhile  to  illustrate  the  synthesis  of  a 
single  cluster  that  reduces  the  distortion  in  the  secondary  pattern  caused 
by  the  lateral  displacement  of  the  feed.  We  refer,  once  again,  to  the 
pattern  shown  in  Fig.  9.  We  also  plot  the  corresponding  "equivalent 
aperture"  distribution  [see  2.4]  whose  Fourier  transform  gives  the 
secondary  pattern  under  the  so-called  "small  angle  approximation" 
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[see  Ref.  3].  We  plot  this  distribution  in  Fig.  10  for  three  different 
azimuthal  directions,  viz.,  <(i'  =0°,  90°  (same  as  -90°)  and  180°.  It  is 
evident  that  the  phase  error  is  quite  large  in  the  (p'  = 0°  and  180° 
directions.  We  point  out  that  the  error  in  the  phase  is  the  departure  from 
the  linear  taper  which  only  shifts  the  direction  of  the  beam.  The  cluster 
excitation  coefficients  that  minimize  the  phase  error  are  found  from  (4.8) 
and  (4.19)  and  are  given  by 


and 


(7.4676  + j 0.7144)  x 10 


(9.8688  + j 1.2031)  x 10 
(1.0088  + j 0.1170)  X 10 


(0.0709  + j 0.0099)  x 10"^  J 


-0.1032  - j 0.0147 
Ij  = -0.0334  - j 0.0060  I 
Ig  = 0.1626  + j 0.0179  I 
0.2929  + j 0.0336  , 


(4.29) 


(4.30) 


Fig.  10.  Equivalent  aperture  dis- 
tribution for  a single 
feed  who,e  center  ele- 
ment is  displaced  3x  in 
the  -X  direction. 


It  is  interesting  to  compare  the  com- 
pensated and  uncompensated  phase 
distribution  in  order  to  assess  the 
degree  of  improvement  achieved  by  using 
a cluster  feed  as  opposed  to  a single 
one. 

We  have  plotted  the  compensated 
case  in  Fig.  11  which  should  be  com- 
pared with  the  original  phase  distri- 
bution shown  in  Fig.  10.  The  following 
features  may  be  noted  from  Fig.  11. 

The  amplitude  taper  has  increased  from 
-10  dB  to  -14  dB  in  the  ([i'  = 0 and 
180°  plane  and  to  -12  dB  in  the 
p'  = 90°  plane.  The  improvement  in 
the  phase  error  is  notable  in  the 
range  0 < s < 0.8,  where  the  amplitude 
distribution  is  relatively  significant. 
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Fig.  11.  Equivalent  aperture  dis- 
tribution for  a cluster 
feed  whose  center  ele- 
ment is  displaced  3x  in 
the  -X  direction. 


In  Fig.  12  we  plot  the  compensated 
secondary  pattern  which  should  be  com- 
pared to  the  corresponding  uncompensa- 
ted pattern  in  Fig.  8.  Although  the 
compensated  pattern  does  not  by  any 
means  reproduce  the  on-focus  feed 
pattern,  it  nevertheless  reduces  the 
distortion  in  the  secondary  pattern 
quite  noticeably. 

The  sidelobe  level  in  the  a =-180° 
direction  is  reduced  from  -13  dB  to 
-17  dB  and  the  shape  of  the  pattern 
improves  in  the  a = 0°  plane. 

Numerical  experiments  have  been 
carried  out  by  varying  the  parameter 
and  the  weight  function  W(s)  of 
the  minimization  integral;  the  basic 
findings  are  that  the  method  works  as 


(«)  (b)  (c) 


Fig.  12.  Compensated  secondary  pattern  with  a cluster  feed  synthesized 
using  equivalent  aperture  method  analytical  solution. 
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i well  as  with  an  arbitrary  and  that  varying  the  weight  function  has 

relatively  little  effect  on  the  results. 

Having  derived  the  cluster  coefficients,  we  can  use  them  for  synthe- 
sizing contour  beams  in  a relatively  straightforward  manner  if  we 
‘ choose  to  use  the  spotlighting  method  for  synthesizing  the  contour 

coverage.  In  this  method,  one  merely  uses  a sampling  of  the  specified 
contour  pattern  in  the  secondary  to  determine  the  excitation  coefficients 
of  the  individual  cluster.  Note  that  in  this  procedure  an  implicit 
assumption  is  made  on  the  behavior  of  the  phase  of  the  secondary  pattern 
which  is  not  specified.  Consequently,  we  do  not  expect  the  results  thus 
derived  to  necessarily  be  optimum  in  the  sense  of  a "best-fit"  criterion 
imposed  on  the  specified  power  pattern.  However,  the  cluster  spotlighting 
method  is  not  only  extremely  efficient,  but  the  result  obtained  via  the 
use  of  this  approach  often  are  close  to  optimum.  For  this  reason  the 
results  derived  via  the  spotlighting  technique  can  also  serve  as  good 
starting  points  for  other  iterative  techniques. 

In  Fig.  13  we  show  the  ETZ  contour  as  viewed  from  a synchronous 

orbiting  satellite  located  at  the 
equator  and  at  the  longitude  of 
112°  W.  The  scale  for  this  figure 
is  in  kau  and  kav  with  D/X  = a/2x 
= 151.05.  A reflector  of  this  size 
produces  a beam  of  -3  dB  beamwidth 
as  depicted  in  the  inset  in  Fig.  13 
when  the  primary  feed  is  described  by 
(4.8)  with  n - 2.254  that  gives  a 
-10  dB  edge  taper.  Recall  that  an 
on-focus  feed  illuminating  this 
reflector  generates  a secondary 
pattern  with  a -26  dB  sidelobe. 

An  equilateral  triangular  lattice 
of  beams  is  superimposed  in  Fig.  13 
showing  the  locations  of  the  centers 
of  the  spotlighting  beams.  The 
separation  distance  between  the  feed 
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Fig.  13.  Eastern  Time  Zone  - 
synchronous  satellite 
view  from  equator  at 
1120  W longitude. 

0 = 151. 05X. 


elements,  which  Is  a design  parameter  for  the  system,  is  chosen  to  be  1, 
which  produces  a -6  dB  intersection  level  between  adjacent  beams.  Such  an 
arrangement  assures  minimum  ripple  magnitude  when  well-formed  beams  are 
added  in  space. 

For  the  size  of  the  reflector  chosen  for  the  example,  only  nine 
clusters  can  be  positioned  such  that  their  main  peaks  fall  inside  the 
contour  which  spans  about  2°  as  viewed  from  the  synchronous  orbit.  A 
total  of  26  beams  are  associated  with  these  nine  clusters. 

The  resulting  contour  beam  realized  with  the  nine-cluster  feeds  is 
shown  in  Fig.  15  and  the  desired  power  pattern  of  the  contour  is  plotted 
in  Fig.  14.  The  entire  design  can  be  carried  out  using  only  a few 
seconds  of  computer  time  when  the  cluster  superposition  method  is  used. 

The  field  calculation,  even  at  nearly  ten  thousand  or  so  observation 
points,  can  be  carried  out  in  a rapid  manner  using  the  series  approach 
described  in  Section  2.  It  is  evident  from  Fig.  15  that  the  resolution 
achieved  from  the  use  of  a 150  A diameter  dish  is  not  sufficient  to 
accurately  follow  the  contour.  Increased  resolution  can  be  obtained  by 


Fig.  14.  ETZ  objective  power  pattern 
function  |Eg(u,v)|2. 

ka  = ttO/A  = IT  15 1.05. 


Fig,  15.  ETZ  contour  power^pattern 
function  lE(u,v)p. 
ka  * ttD/a  * wl51.05. 


i 

t 

t 


I 


using  a larger  dish;  we  show  in  Fig.  16  the  results  for  a 302.1  X dish. 
Note  that  the  diagrams  have  been  scaled  by  the  radius  a,  which  is  now 
twice  as  large  compared  to  the  one  in  the  last  example.  The  number  of 
clusters  that  are  accommodated  within  the  contour  now  increases  to  32  and 
the  associated  number  of  single  beams  is  63.  The  resolution  of  the 
synthesized  contour  pattern  is  substantially  improved.  It  should  be 
pointed  out,  however,  that  beam  distortion  puts  a severe  limit  on  the 
maximum  achievable  resolution  and  one  cannot  expect  continuous  improvement 
with  an  increase  of  reflector  size.  Although  not  apparent  at  first  sight, 
the  increase  in  computation  time  due  to  the  increase  in  the  number  of 
beams  is  very  little.  This  is  a significant  advantage  of  the  cluster 
synthesis  approach  developed  in  this  section. 

4.4  Gradient  Optimization  Method 

In  the  last  section,  we  discussed  a procedure  for  synthesizing  a 
contour  beam  using  a spotlighting  technique  to  determine  the  excitation 
coefficients  of  a set  of  cluster  feeds.  The  tactic  employed  in  this 
approach  was  to  closely  approximate  an  ideal  "equivalent  aperture  distri- 
bution" with  that  produced  by  each  cluster  by  varying  the  relative 
amplitudes  and  phases  of  the  individual  elements  of  the  cluster.  In 
following  this  procedure,  an  implicit  assumption  is  made  that  a close 
approximation  to  the  ideal  distribution  in  the  "equivalent  aperture" 
will  produce  a relatively  distortionless  secondary  (field)  pattern. 
However,  no  convenient  quantitative  correlation  between  the  "goodness" 
of  the  distribution  in  the  equivalent  aperture  and  the  lack  of  distortion 
in  the  secondary  power  pattern  can  be  readily  established.  One  can, 
therefore,  choose  to  work  directly  with  the  secondary  power  pattern  and 
seek  to  approximate  the  same  using  some  type  of  minimization  procedure. 

To  this  end  one  can  define  a performance  index  A given  by 


f'2 

^2 

9 

A = 

du  1 

dv  W(u,v) 

|E(u,v)|  - |EQ(u,v)r 

V 

^1 

w w 

where  W(u,v)  is  a weighting  function,  E(u,v)  is  calculated  from  (4.1)  and 
Eq  is  the  desired  (real)  power  pattern  which  can  be  taken  to  be  the 
undfstorted  pattern  for  a symmetric  feed.  The  region  of  integration  in 
(4.31)  should  be  taken  to  be  an  area  in  the  (u,v)  space  which  is  large 
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enough  to  cover  the  range  of  interest  in  the  secondary  pattern  space.  Any 
number  of  available  optimization  schemes  can  be  employed  to  minimize  A and 
derive  the  solution  for  the  excitation  coefficients  I^'s  that  implicitly 
appear  in  (4.31)  through  the  expression  for  E(u,v).  The  numerical  results 
for  Ij^'s  obtained  from  an  optimal  solution  of  (4.31)  are  given  by 


I2  = -0.0816  - j 0.0734 
12=17=  “O-O^SO  + j 0.0141 
= ig  = 0.2034  - j 0.0290 
Ir  = 0.2833  - j 0.1771 


(4.32) 


which  should  be  compared  with  the  one  derived  by  working  with  the 
"equivalent  aperture"  method  and  given  in  (4.30).  The  secondary  pattern 
for  this  choice  for  I^'s  is  plotted  in  Fig.  17.  We  note  that  at  least  for 
this  example  the  optimal  solution  is  not  significantly  superior  to  the 
secondary  pattern  derived  earlier  using  the  "equivalent  aperture"  method. 
However,  the  computer  time  required  for  deriving  the  optimal  solution  is 
considerably  higher  than  that  necessary  to  extract  the  cluster  coefficients 
using  the  procedure  outlined  in  Section  4.3.  We  also  mention  that  the 
solution  for  derived  from  the  "equivalent  aperture"  approach  can  serve 
as  an  excelleit  initial  guess  in  the  gradient  optimization  program. 

Finally,  it  is  possible  to  employ  the  optimization  approach  to  synthesize 
the  entire  contour  beam  directly,  rather  than  just  the  constitutent  cluster 
beams.  It  is  not  difficult  to  see,  however,  that  the  time  involved  in 
deriving  an  optimal  solution  directly  for  the  entire  contour  can  be  very 
large  and  that  there  is  the  increased  risk  of  converging  to  a false  minimum. 

Although  we  have  not  discussed  the  algorithmic  details  of  optimization 
procedures,  we  mention  that  a number  of  optimization  subroutines  are  usu- 
ally available  in  a typical  computer  program  library.  More  sophisticated 
procedures  for  minimization  techniques,  such  as  the  minimax  procedure, 
may  be  found  in  [9].  An  application  of  this  method  to  the  contour  beam 
synthesis  problem  has  been  discussed  in  [10]. 
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Appendix  I. 


Beam  Deviation  Factors  for 
Symmetric  Reflectors  (Parabolic) 

The  beam  deviation  factors  B , B are  defined  by 

u V 


B 

u 


B 

V V3 


(A. 1.1) 


where  we  recall  that  and  Vj^^  are  the  direction  cosines  of  the  beam 
maximum  and  and  v^  are  given  by  = e^/f  and  v^  = e^/f.  The  beam 
deviation  factors  are  generally  independent  [11]  of  the  feed  displacement, 
at  least  when  and  v^  are  not  very  large.  It  is  also  fortuitous  that 
these  parameters,  which  are  essentially  constants,  can  be  found  directly 
from  the  primary  pattern  of  the  feed  centered  at  the  focal  point  by  using 


^ (1  - cos  0')[3Hp  + Ep]  ^ ds 


u rl 


(1  - cos  0')[3Hp  + Ep]  s-^  ds 


(A. 1.2) 


where  Ep  and  Hp  are  the  E-plane  and  H-plane  patterns  of  the  feed.  The 

other  constant  B^  can  also  be  derived  from  (A. 1.2)  by  simply  interchanging 

E^  with 
P P 


Appendix  II. 


A Method  for  Computing 
the  "Equivalent  Aperture"  Distribution 


We  will  now  briefly  outline  a method  for  computing  the  "equivalent 
aperture"  distribution  for  a given  feed  whose  primary  pattern  is  known. 
Let  the  feed  pattern  be  specified  in  the  (e',((>')  coordinate  system 
(Fig.  1)  as 


t - ® 

^feed 


-jkp 


— t®'  Eg.  E..] 


(A. 2.1) 


Assuming  that  the  E-  and  H-plane  patterns  of  the  feed  can  be  smoothly 
interpolated,  we  obtain  for  a nominally  y-polarized  feed 


E , H E (e')  sin  <p'  , E , = -H  (e')  cos 

up  ((>  p 


(A. 2. 2) 


so  that 


^('^.'1'')  = - (1  - cos  e')  (-Hp  + Ep)  sin  (t>'  cos  4>' 


} H 

[_l 


+ y|Hp  cos  (^'  + Ep  si 


in^  ij)' 


+ z 


e' 


ctn  -y  sin  (p'  E 


(A. 2. 3) 


Each  rectangular  component  of  ? can  now  be  related  to  a corresponding 
component  of  the  surface  current. 
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